Mind the gap but also the spin: why the Heyd-Scuseria-Ernzerhof hybrid functional 

description of VO2 phases is not correct 
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In contrast with recent claims that the Heyd-Scuseria-Ernzerhof (HSE) screened hybrid functional 
can provide a good description of the electronic and magnetic structure of VO2 phases [V. Eyert, 
Phys. Rev. Lett. 107, 016401 (2011)], we show here that the HSE lowest-energy solutions for both 
the low-temperature monoclinic (Ml) phase and the high-temperature rutile (R) phase, which are 
obtained upon inclusion of spin polarization, are at odds with experimental observations. For the 
Ml phase the groundstate is (but should not be) magnetic, while the groundstate of the R phase, 
which is also spin-polarized, is not (but should be) metallic. The energy difference between the 
low-temperature and high-temperature phases is also in strong discrepancy with the experimental 
latent heat. 



The screened hybrid functional approach of Heyd, 
Scuseria and Ernzerhof (HSE) has accumulated signif- 
icant success in the description of structural and elec- 
tronic properties of molecules and solids at a moderate 
computational cost PHI]- It has been recently argued 
that density functional theory (DFT) calculations based 
on this functional are well capable of describing the elec- 
tronic and magnetic properties of the metallic and insu- 
lating phases of vanadium dioxide VO2 [4]. In particular, 
it was shown there that HSE calculations are capable of 
producing the expected band gap in the electronic struc- 
ture of the monoclinic phase, in contrast with other DFT 
approximations [5j]. This is exciting progress, because 
modelling the VO2 phase transition with relatively inex- 
pensive DFT methods (the alternative is many-body GW 
calculations or dynamic mean- field theory |6j-l8j) opens 
the door to a more active role of ab initio design in 
the development of applications such as "smart" ther- 
mochromic windows [9j. 

When heated to ^340 K, pure VO2 exhibits a tran- 
sition from a monoclinic (Ml) semiconductor phase to 
a tetragonal, rutile- like, metallic phase (R) [TOj , The 
transition is first-order [ll| , and a latent heat of 44 meV 
per VO2 formula unit has been measured by calorimet- 
ric methods [12]. Therefore it can be expected that ac- 
curate calculations yield a band gap for the Ml phase 
but not for the R phase, and that the calculated total 
energy per formula unit of the R phase is higher than 
that of the Ml phase. In order to compare the electronic 
structures and energies of the two phases we have per- 
formed HSE calculations using the planewave DFT pro- 
gram VASP [lH, first using non-spin polarized calcula- 
tions as in Q , and then using spin-polarized calculations 
in different magnetic configurations. Experimentally de- 
termined crystal structures were used in the calculations 
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14], including four and two formula units for the Ml 
and R phases, respectively, without geometry relaxation. 
In order to ensure a reliable energy comparison between 
phases, precision parameters were chosen to achieve a 
convergence of 1 meV/atom in total energy. This re- 
quired k-point grids of 5x5x5 and 6x6x9 for the Ml and 
R unit cells, respectively (the same grids were used for the 
non-local exchange contributions). The energy cutoff for 
the planewave basis set was 400 eV, and the interaction 
between the core (up to 3p for V and up to Is for O) and 
valence electrons was described with the projected aug- 
mented wave (PAW) method [HI . The standard settings 
were used for the screened hybrid functional, i.e. 25% 
of Hartree-Fock exchange was mixed in, with a screen- 
ing parameter fx — 0.207A -1 , and the local contributions 
were calculated with the Perdew-Burke-Ernzerhof (PBE) 
functional [16j . 

The results for the non-magnetic (NM) calculations in 
Fig. [1] show that the energy of the R phase is indeed 
higher than that of the Ml phase. However, the calcu- 
lated energy difference (£[R]-£'[M1]=232 meV per V0 2 
formula unit) between the two phases is too large com- 
pared with the experimental latent heat. Our estimation 
of the latent heat ignores the phonon contribution, but 
this can be expected to be relatively small (the difference 
between the zero point energies of the R and Ml phases 
has been estimated to be -17 meV per formula unit using 
shell model calculations 17]). Furthermore, upon inclu- 
sion of spin polarization in the calculations, the total en- 
ergies for the R and Ml phases are lowered with respect 
to the respective NM solutions, which implies that the 
HSE electronic groundstates for both phases are not the 
ones described in [||. The agreement between HSE-level 
theory and experiment should then be revised. 

Allowing spin polarization in the Ml phase calculation, 
with either antiferromagnetic (AFM, with magnetic mo- 
ments alternating orientations along the rutile c axis) or 
ferromagnetic (FM) configurations, lowers the total en- 
ergy by 463 meV and 365 meV per formula unit, respec- 
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TABLE I. Band gaps obtained from HSE calculations of VO2 
phases in different magnetic configurations. Experimental 
values are from 0,1. 
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FIG. 1. Energies per VO2 formula unit of the monoclinic 
(Ml) and rutile (R) phases of VO2 as obtained with the HSE 
functional in nonmagnetic (NM) , ferromagnetic (FM) and an- 
tiferromagnetic (AFM) configurations. All energies are given 
with respect to the NM solution for the Ml phase, which 
should be (but is not for HSE) the global groundstate at zero 
temperature. The experimental latent heat is taken from ref- 
erence Hal. 



tively, in comparison with the NM calculation. Magnetic 
moments of ~1 \ib arc found by integrating the spin den- 
sity around the V ions, both in the AFM and in the FM 
case. As shown in Table U and illustrated in the density of 
states (DOS) plots of Fig. [3J these lower-energy magnetic 
solutions have wider band gaps than the NM solution, 
thus worsening the agreement with experiment. More 
importantly, the existence of a groundstate with local 
magnetic moments is in conflict with the well-established 
non-magnetic character of the Ml phase. Experimen- 
tally, VC>2(M1) exhibits only a very small, positive and 
temperature-independent magnetic susceptibility, which 
is associated with van Vleck paramagnetism, and there- 
fore has no magnetic moments in the groundstate (l8l . [l9| . 
The wrong groundstate appears in HSE calculations due 
to the Hartree-Fock exchange mixed in the HSE func- 
tional, which is required for opening the band gap in 
the Ml phase [J], but at the same time tends to stabi- 
lize localized magnetic moments, in this case excessively 
(although in other systems the same effect can actually 
lead to better agreement with experiment [20|). In our 
calculations we are using the experimental crystal struc- 
ture, but we have checked that the magnetic moments 



on the V ions remain stable in the HSE solution for the 
Ml phase upon changes of +/- 1% in the cell parame- 
ters (HSE-optimised lattice parameters typically deviate 
less than 1% from experimental values J3[). We find that 
the difference between the non-magnetic and the mag- 
netic solution (taken with FM ordering for calculation 
convenience) increases 13% when the cell is expanded, 
and decreases 13% when the cell is compressed, but the 
solutions with magnetic moments are always more stable 
than the non-magnetic ones. 

For the R phase, HSE gives a ferromagnetic (FM) insu- 
lator groundstate (Fig. [2]), with local magnetic moments 
of ~1 [Ib on the V atoms, which is 739 meV per formula 
unit below the NM solution. This solution is also 506 
meV lower in energy than the NM solution for the Ml 
phase. The antiferromagnetic (AFM) solution, with V 
magnetic moments in alternate orientations along the c 
axis, also has a band gap, and is close in energy to (but 
less stable than) the FM groundstate. This situation is 
again clearly unsatisfactory. Not only it is wrong from 
a thermodynamic point of view that the groundstate for 
the R phase is significantly lower in energy than the Ml 
phase, but also the magnetic insulator groundstate found 
by HSE for the VC>2(R) structure is in conflict with ex- 
perimental evidence. The problem here is not the mag- 
netic character, as VC>2(R) seems to be paramagnetic 
with a large temperature-dependent susceptibility [l9j . 
but the insulating character of the solution. HSE pre- 
dicts large band gaps for both the FM and AFM solutions 
(Table Q}. This is a serious problem, because the metal- 
licity of the R phase is well-established experimentally 
by conductivity measurements (e.g. [HI, [Hj]). A note of 
caution should be added here: our spin-polarized calcu- 
lations of the R phase are based on magnetically ordered 
cells, while the proximity in energy between the FM and 
AFM solutions, and also the experimental evidence, sug- 
gest that the system should have paramagnetic disorder. 
It is therefore still possible in principle (although rather 
unlikely based on the calculated gap values for the or- 
dered configurations) that the HSE can recover the cor- 
rect metallic solution for the R phase if magnetic disorder 
could be accounted for, due to broadening of the bands. 
Unfortunately, this type of calculation at the HSE level 
is beyond our computing capabilities at the moment. 

We also note that the HSE approximation to the latent 
heat is pretty bad, regardless of which solutions are taken 
for the calculation. The difference in energy between the 
magnetic R phase and the non-magnetic Ml phase is very 
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FIG. 2. Electronic density of states (DOS) of the nonmagnetic (NM), ferromagnetic (FM) and antiferromagnetic (AFM) HSE 
solutions for the monoclinic (Ml) and rutile (R) phases of VO2. 



large and of opposite sign compared to experiment. If 
we take the (physically wrong) ground-state solutions for 
both phases, the absolute value of the HSE latent heat 
becomes much lower, but still has the wrong sign. This 
is interesting, because previous calculations based on the 
local density approximation (LDA), which is unable to 
account for the band gap opening in VOa(Ml), did show 
good agreement (within 10 meV) with experiment in the 
latent heat of the transition [2J|. The HSE functional 
thus performs worse than the simple LDA functional in 
the description of the relative energies of the phases. 

The existence of incorrect groundstates severely limits 
the usefulness of the HSE approach in the investigation 
of this important oxide and its phase transitions. Al- 
though it can be argued that in the case of the Ml phase 
one can still obtain a meaningful solution by forcing a 
non-spin-polarized calculation, such an approach is not 
satisfactory. To illustrate why, we consider the case of the 
tungsten-doped VC>2(M1) phase, which is interesting for 
applications (2 at.% doping with tungsten can lower the 
scmiconductor-to-metal transition point to room temper- 
ature i24j). We have investigated the doped oxide using 
a 96-atom supercell with one V substituted by W. Since 
the substitution leads to an odd number of electrons per 
cell, a non-spin polarized calculation necessarily leads to 
a metallic solution, and in fact this will be the case re- 



gardless of the W concentration employed in the simula- 
tion. This is an artifact of the spin-restricted calculation; 
experimental observations confirm that W-doped VO2 re- 
mains a semiconductor at low temperatures [2~il |25| . In 
order to deal correctly with the odd number of electrons, 
spin polarized calculations are necessary; however, they 
lead to localized magnetic moments in all the V ions in 
the cell, i.e. to the wrong groundstate. This was the 
case even when the initial magnetic moments of the V 
ions were set to zero in the calculations. Therefore, the 
problem of W doping in VC>2(M1) becomes intractable 
within the HSE approximation. 

In summary, although the HSE description of the band 
gap opening in V02(M1) reported recently [4] is wel- 
comed progress, the results presented here show that the 
HSE functional does fail in the description of both the 
electronic structure and the energetics of the transition of 
VO2. It gives a magnetic groundstate for the Ml phase, 
a non-metallic groundstate for the R phase, and an R-Ml 
energy difference in significant disagreement with the ex- 
perimental latent heat. Vanadium dioxide thus continues 
to be a challenge to band theory. Despite its succesful 
record, the HSE functional needs to be used cautiously, 
particularly in the simulation of the magnetic properties 
of transition metal oxides. 
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